On Slow Light as a Black Hole Analogue 
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Although slow light (electromagnetically induced transparency) would seem an ideal medium in 
which to institute a "dumb hole" (black hole analog), it suffers from a number of problems. We show 
that the high phase velocity in the slow light regime ensures that the system cannot be used as an 
analog displaying Hawking radiation. Even though an appropriately designed slow-light set-up may 
simulate classical features of black holes - such as horizon, mode mixing, Bogoliubov coefficients, 
etc. - it does not reproduce the related quantum effects. 
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I. INTRODUCTION 

The astonishing ability to slow light to speeds of a 
few meters per second has been a striking development 
in quantum optics, see e.g. [1]. The idea to use matter 
systems as analogs [2] to the (yet unobserved) Hawking 
effect [3] for black holes has raised the possibility of ex- 
perimentally testing certain assumptions which enter into 
those calculations, see e.g. [4]. The dependence of those 
analogs on the detection of sound waves however causes 
problems, as the detection technology for light is much 
more developed than for sound, and finding an optical 
analog to black holes [5-8] could make the experimen- 
tal detection of the analog for Hawking radiation easier, 
cf. [9]. 

Recently Leonhardt [5,7] has suggested that slow light 
systems could be used to create such an analog, but that 
approach has been criticized by one of us [9] . This paper 
is an amplification of that criticism, looking in detail at 
the use of slow light in such an analog, and trying to un- 
derstand in what sense slow light could be used to create 
and analog for black holes, and why, despite that analog, 
it will not create the thermal radiation characteristic of 
the Hawking process. 



II. DESCRIPTION OF THE SET-UP 

In order to generate slow light, one first chooses 
an atom with a convenient set of atomic transitions, 
cf. [1,10]. In particular, a system is chosen with two 
long lived meta-stable or stable states, and with one 
state which is coupled to these two states via dipole elec- 
tromagnetic transitions (A-system). Let us call the two 
lower meta-stable states \a) and |6). The third higher en- 
ergy state is |c). The two states \a) and \b) are assumed 
to have energy — oj a , —ujh, and |c) has energy zero and 
decay constant L > 0. (I.e., this higher energy state is 
assumed to have decay channels other than electromag- 
netic radiation to the |a) and \b) states.) 

The electromagnetic field, which we will assume has 
a fixed polarization, will be represented by the vector 



potential A where E = dtA (temporal gauge). 



A. Effective Lagrangian 

The effective Lagrangian for this system can be written 
as (H = c = 1 throughout) 



L 



dx C J 



E(4 



-Aip 



(1) 



with the usual term governing the dynamics of the elec- 
tromagnetic field 
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C A = ~ [E 2 



1 



B 2 ] = \ [(d t Af - {d x Af] 
and the Lagrangian of the atomic states 



(2) 



L j = i (V'aj & 1paj + ?Pb3 9 t Tpbj + 1p*j 9t i>, 



(3) 



as well as the interaction term in dipole approximation 



-Aiji _ 



E(xj) (e a ipejipaj + £fc 4>* c jVJbj) + h 



(4) 



where Xj is the location of the j-th atom. Here the 
are the amplitudes for the j-th particle being in the cor- 
responding state |atom j) — ip a j \a} + y^ b j \b) + ij) c j |c) and 
e a , £{, are the associated dipole transition amplitudes. 

In contrast to the usual set-up, i.e., a strong con- 
trol beam and a weak (perpendicular) probe beam, let 
us assume that there is a strong background counter- 
propagating electromagnetic field 



A (t,x) 



sin f 



a iu>b(t+x) 



h.C. 



(5) 



i.e., at the resonant frequencies of the two transitions. 
The mixing angle 8 controls the relative strength of the 
left- and right-moving beam and Q denotes the averaged 
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Rabi frequency*. For a single beam (9 = or 9 — 7r/2) Q, 
reduces to the exact Rabi frequency of that beam. The 
fact that the phase velocity is unity (i.e., the light speed) 
prefigures the fact that the effective dielectric constant of 
the atoms is unity at these transition frequencies when 
the atoms are in the so called "dark state", cf. [1,10]. 

In the following we shall assume that we can and arc 
making the rotating wave approximation. One solution, 
the only (up to an overall phase) non-decaying solution, 
for the atoms is 



raj 
<j 

€j 



-°. = -e 

V& = o. 



,ia,.(t-x J ) si]lfl} 

cos 0, 



(6) 



Since the Rabi oscillations between the states \a) and 
|c) interfere destructively with those between the states 
\b) and |c) (leading to a vanishing occupation of |c)), this 
solutions is called a dark state (no spontaneous emission) . 



B. Linearization 

Let us redefine our electromagnetic field such that 
A(t, x )=(n — + <i> a (t, x)) e -*".(*-*) 



sin 9 



eb^b 



+ Mt, *)) e-^ t+x) +h.c, (7) 



where we are going to assume that both Q b and $ a are 
slowly varying functions of time and space (i.e., beat fluc- 
tuations). 

Furthermore, let us define 



= (*»,•+ sin 0) e^^'), 
^bj = (* tj --cos0) e^* <«+*'>, 
i>cj = *cj , 



(8) 



where the new variables \& are also assumed to be slowly 
varying. 

Substituting into the Lagrangian, retaining only the 
second order terms''' in the 'J, $b, $ a , using the rotating 
wave approximation, and neglecting time derivatives of 
$b and $ with respect to uj a and u> b we get the effective 
(approximated) Lagrangian for the beat fluctuations 

C A ~ 2icu a <t>* a (8 t + d x )$ a + 2iuo b ^l(d t - d x )$ b , (9) 



*Note that is often defined differently, i.e., with an addi- 
tional factor of two. 

^The zeroth-order contributions decouple and the first-order 
terms vanish after an integration by parts, since the back- 
ground fields solve the equations of motion. 



and the atomic states 

if ~ i (*v d t * aj + ** bj d t * bj + K 3 d t * cj + r K^c 3 ) 

-iSl (V^Vaj cos 9 + ** cj * b j sin 9 - h.c.) , (10) 
as well as the interaction 

Lf 9 ~ -iw a e a sme^ a (xj)^* cj + iw b e b cos 9 <f> b (xj)^* cj 



+h.c. 



(11) 



III. EQUATIONS OF MOTION 

The equations of motion for the particle amplitudes 
can be derived from the effective Lagrangian 

d t ^ aj = -Q cos9^ cj , 
dt^bj = -ft sin0tf cj , 
dt^cj = ft (cos 9 + sin 9 * w ) - T 

+w a e a sin0$ o (a;j) - W{,e fc cos0<&6(a;j) , (12) 

and the equation of motion for the fields <& a and $b are 

2 (St + 9 X ) $ a = -e a sin * CJ - 5 (a; - Xj) , 

j 

2 (d t - 9,) $ 6 = +e fc cos J2 *d 5 ( x ~ x o) ■ (I 3 ) 

j 

Assuming that the particles are sufficiently closely spaced 
so that there are many particles in a space of the order of 
a wavelength of the field, the sum over j can be replaced 
by the density of the particles 

2 (d t +d x )$ a = -p(x) e a sin 9 t> c (x) , 
2(d t -d x )$ b = +p{x)e b cos0* c (x). (14) 



A. Effective Dispersion Relation 

Assuming harmonic space-time dependence e -" lut + lKX 
of all of the variables, we can solve the equations of mo- 
tion for the atomic amplitudes (12) 

*cj( W ) = \ e a^a SOI 6 <& a (uJ,Xj) - e b Ul b COS 6 Q b {u), Xj)] 



LU 2 -Q 2 +iTLU 



(15) 



and inserting this result into Eq. (14) we finally obtain 
the dispersion relation 

(lu + X{u) -k)(uj + Y(lu) +k)= X{uj)Y(lu) , (16) 

where 

u puj a e 2 a sin 2 9 

- 77 



2 n 2 - oj 2 - iTuj ' 
_ w pw b el cos 2 9 
[UJ) _ 2 tt 2 -uj 2 -iTuj 



(17) 
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B. Adiabatic Regime 

For small u and k, the dispersion relation derived 
above turns out to be linear, i.e., lo oc k. Let us spec- 
ify the required conditions. As already mentioned above, 
Eq. (14) is valid for wavelengths which are much larger 
than the inter-atomic distance Ax (typically a few hun- 
dreds of nanometers) only 



k <C 



Ax 



(18) 



In addition, the manipulations of the previous Section 
(rotating wave approximation) are based on the assump- 
tion that the fields $b and $ a are slowly varying, i.e., 
u <C u) a , &b- However, since the Rabi frequency £1 is 
supposed to be much smaller than the atomic transition 
energies u) a , u) b and the decay rate is assumed to be small 
r < 0, the knee frequency ft of the above dispersion re- 
lation yields the relevant frequency cut-off 



f ^ 2 1 

u> < min < Q, uj a , w , — > = Q . 



(19) 



In this limit, i.e., in the adiabatic regime, Eq. (12) can 
be solved via 



io a £a sin V ■ u b e b cos V ■ 

W n = — $ a ®b ■ 



f! 2 

Rescaling the fields via 



ft 2 



$a = w a e a sin #$ a , 
$ b = uj b e b cos 8 $ b , 



(20) 



(21) 



Eqs. (13) and (14) become 

pu a el sin 2 6 



(d t + d x )$ a = - 
(d t - d x ) $ b = + 



2n 2 

p LO b e 2 cos2 
202 



(dt$ a - dt$ b ) 

(c\$ - d t $ b 



(22) 



In order to cast these two first-order differential equations 
into the usual second-order form, let us choose 8 such 
that* 



puj a e 2 sin 2 8 p^ b ^ 2 cos 2 8 



20 2 



20 2 



= H . 



(23) 



where the dimensionless quantity H describes the slow- 
down of the waves and can be very large H 3> 1. In 
terms of the fields 



(24) 



* Otherwise one would obtain an velocity- like term even for 
a medium at rest, cf. Sec. V below. However, this term alone 
cannot generate an effective horizon. 



we can indeed combine the two first-order equalities 
above into one second-order equation 



dt 2 



i) 



1 



d 



1 + 2K dx 



$+ = 0. 



(25) 



Obviously, small background fields, i.e., small Rabi fre- 
quencies f2, may generate a drastic slow-down N> 1. 

Note, however, that the above wave equation differs 
from the equation of motion describing a slow-light pulse 
in the usual set-up - i.e., a strong control beam and a 
weak (perpendicular) probe beam, cf. [1,10] 



([i + N]a t ±a x )* = o. 

Hence the slow-down in Eq. (25) 



(26) 



^group 



Wpha 



l/Vl + 2H of the design proposed in the present arti- 
cle is not as extreme as that of the usual set-up i> gr0 up = 
1/(1 + N) 7^ w p hasc ~ 1, but still substantial. 



IV. EFFECTIVE GEOMETRY 

So far we considered a static medium at rest with a 
possibly position-dependent H = H(x). Now we allow 
for a space-time varying variable H = N(t, x), where the 
medium is still at rest. A change of H can be generated 
by varying p, i.e., by adiabatically adding or removing 
atoms. The other parameters in Eq. (23) remain con- 
stant - a time-dependent £1, for example, would generate 
additional source terms and thereby invalidate the back- 
ground solution. 

Furthermore, we shall assume u a — u b as well as e a — 
e a (which is a reasonable approximation) and hence 8 = 
7r/4 for the sake of simplicity and absorb these quantities 
by rescaling the fields $±. 



A. Effective Action 

Introducing the abbreviation * = (^a, <!/{,, , J C ) T the 
linearized Lagrangian governing the dynamics of the *5f- 
fields in Eqs. (10) and (11) can be cast into the following 
form 



A 9 = J d 2 x (W f * + * f M * 



+ 



■)■ 



(27) 



with M denoting a (self-adjoint) 3x3 matrix and N a 
three-component vector as determined by Eqs. (10) and 
(11). In terms of the differential operator defined via 

D = idt + M and its formal inverse D we may com- 
plete the square 

A 9 = J d 2 x (\J> f • D ■ * - <P*_ ATt • D 1 • AT<I>_) , (28) 



3 



with 



* = * + • AT$_ 



(29) 



Assuming that the quantum state of the ^-fields is ad- 
equately described by the path-integral with the usual 
(regular) measure we are now able to integrate out 
(i.e., eliminate) those degrees of freedom explicitly arriv- 
ing at an effective action for the $-fields alone 



exp{i4 eff } = --^ / 3* cxp{z(^*+^*)} 



(30) 



As demonstrated in Eq. (28), the above path-integral is 
Gaussian (ID* = 5)*) and the associated Jacobi deter- 
minant is independent of <&. Hence we obtain 



A eS = A* 



d 2 x 



<E>* iV f • D 1 • 



(31) 



As usual, the inverse differential operator D causes 
the effective action to be non-local (in time) - but in the 
adiabatic limit w -C fi, Ho> <C w ,&j K <^ 1/Ax, and <C 
uj a ,b the low-energy effective action is local An 
easy way to reproduce this result is to remember the 
original equation of motion 



D * + iV$ = * = -D -7V$_ 



(32) 



and its solution in the adiabatic limit as given by 
Eq. (20). Together with Eq. (9) we finally arrive at 



+h.c. 



(33) 



Strictly speaking, one obtains an effective action for each 
atom 



A J eS oc i 



dt 



+ h.c. . 



(34) 



where the total effective action incorporates the sum over 
all atoms. With the assumption that the atoms are suffi- 
ciently closely spaced, cf. Eq. (18), and moving in a direc- 
tion perpendicular to the beam (e.g., in the y-direction) 
only, we recover Eq. (33). 

An alternative method for effectively changing the den- 
sity p is to cause transitions between the states \a) and 
| b) and further states \d) and |e), which do not couple to 
the electromagnetic field $± under consideration. The 
dynamics of these additional states is governed by the 
Lagrangian 

£add = #d^d + i4>t i>e + Udlpd^d + Uj e 1pl^e 

+ (i&Ma + K^b + h.c.) , (35) 

where iQ denotes the space-time dependent transition 
amplitude. (This particular parameterization will be 
more convenient later on.) 



If the amplitude (population) of the states \d) and |e) 
is large ipd,e ^ i>a,b and the transition weak Q -C 1, we 
may neglect the back-reaction (ip a ,b — » i>d,e) as wen as 
the associated (quantum) fluctuations and describe the 
process by a classical external source for tp a ,b- 

Furthermore, assuming uj a = uid and uii, = ui e as well 

as 



^ d = + e *".(*-*) sin0 + O(C 2 ), 
?p e = ~e lu "><> t+ ^ cos9 + 0(( 2 ) , 



(36) 



the background solution in Eq. (6) acquires an overall 
pre-factor 



$ = -C(t,x) e ^ t+ ^ cos6, 
# = 0. 



(37) 



This scale factor £(t, x) enters the subsequent formulas 
and effectively changes the density of the contributing 
atoms p e g = Q 2 p. In particular, the wave equation (22) 
gets modified via 



+ d x $+ + 2tt (d t (C$_) = , 



(38) 



which is exactly the same equation as derived from the 
effective action in Eq. (33) with H — > ( 2 H. 



B. Effective Spinor- Representation 

According to Eq. (33) the total effective action for the 
beat fluctuations $± of the electromagnetic field can be 
written as 

A= l -jd 2 x [(1 + 2N) (§*_dt<$>- - [ft** ]*_) 

+ ($*_9 X $+ - [S x $*_]$+)] . (39) 

Introducing the effective two-component spinor i/> (not 
to be confused with the atomic amplitudes ip a ,b,c) 



Vl + 2H$_ 



(40) 



this action can be rewritten as 
d 2 x 



A 



Vl + 2N tydtijp - dt^ ip) 



2J VT+2H 
+ ^a x d x ip - d x ^ a x i/?j - - ^7^ ^iv y ip 



(41) 



with u x , (j y , cr 2 being the Pauli (spin) matrices obeying 



°l = °% = a * = 1 



But this exactly corresponds to the expression for a 
1 +1 dimensional Dirac field x(} 
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= J d 2 x^i 



if we define the Dirac 7-matrices via 

7° = VT+m<j y , 

7 1 = -i<*z , 

and, accordingly, the Dirac adjoint (tptp,tp'y^tp e 

as well as introduce the effective mass 

f 1 an 
to = — - 



2 1 + 2X dx' 
The effective metric is given by {7^,7^} = 2g^ 
dt 2 



ds 2 = 



1 + 2H 



(42) 



(43) 



(44) 



(45) 



(46) 



and displays the expected slow-down. 

For deriving the identity of Eqs. (39) and (42) we need 
the properties of the spin connection (Fock-Ivaneneko 
coefficient) which enters into the spin derivative (remem- 
ber d^ip) = (V^VOV + V'VjWO 

V,W = 3W> + 1>/> , V l ^ = d^-^r tl , (47) 

and is defined by (V^V^V) = (V^h'Y + Vn^V^) 

^ + 1^7" =[7^] , (48) 

with being the Christoffcl symbol. In our 1+1 di- 
mensional representation, the l.h.s. is a linear combina- 
tion of (T y and a z , cf. Eq. (43), and, therefore, the spin 
connection T M has to be proportional to a x . As a result 
we obtain the relation 

{r„,7*} = 0, (49) 

and thus confirm the identity of Eqs. (39) and (42) 

2i£ m=0 = [V„^] t"^ -^7"V^ 

= [d^-f^-^fd^. (50) 

Finally, if we were to choose (over some finite region, 
since H > 0) 



l + 2X(t,x) =f(t)e 



-4mx 



(51) 



the effective mass m in Eq. (45) would be constant 
(which, however, is not necessary for the introduction 
of an effective geometry) and the analogy to the 1+1 
dimensional massive Dirac field complete 



(i7 M V M - to) ip = . 



(52) 



C. Effective Energy 



The energy-momentum tensor of a Dirac field reads 

= - (ip-y(^ v) ip - [V (M V] i V )$) 



= 9 (07(aA)0 ~ [ d (^] 7^)0) , 



(53) 



where the second equality sign holds in general only 
in 1+1 dimensions (in analogy to the simplifications 
above). For an arbitrarily space-time dependent H, how- 
ever, there is no energy or momentum conservation law 
associated to this tensor. But assuming time-translation 
symmetry as described by the Killing vector £ = d/dt we 
may construct a conserved energy via 



E = 



= dx Tq : 



(54) 



which, for the Dirac field in Eq. (53) and the metric in 
Eq. (46), reads 



E = J dx^ (VV' - </> V) 



(55) 



On the classical level, this quantity is (even in flat 
space-time) not positive definite (as is well-known). For 
quantum fields the situation can be different. Impos- 
ing fermionic (i.e., anti) commutation relations the en- 
ergy operator is - after renormalization of the zero-point 
energy and definition of the vacuum state as the filled 
Dirac sea - indeed non-negative (again in flat space- 
time). However, the fields ip = (y/l + 2N $+) T do 
not obey fermionic but bosonic statistics (as one would 
expect, cf. Sec. VI below) and, therefore, the effective 
energy possesses negative parts. 

This fact is not surprising in the context of the elec- 
tromagnetic field, since one has the huge background 
field with which these perturbations can exchange energy. 
However, since in the laboratory frame, the background 
metric is stationary, the energy is a conserved quantity, 
and the potential instability of the negative energy will 
not be triggered. 



D. Inner Product 

Since the (classical) equation of motion (52) can be 
described by means of an effective metric in Eqs. (46) and 
(60) below, we can introduce a conserved inner product 
for two solutions of the wave equation tp 1 and tp 2 - As 
usual, the inner product of the Dirac field can be derived 
by means of the Noether theorem associated to the global 
[/(l)-symmctry ip e lv ip and reads 



W>il^ 2 ) 



dE^ tp 1 7^ ip 2 = / dx 4-9 i>\ 7° '02 



(56) 
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In contrast to the energy in Eq. (55), this quantity is 
non-negative on the classical level (as well as for quan- 
tum fields with bosonic statistics). If we were to impose 
fermionic commutation relations, the above pseudo-norm 
would equal the difference of the number of particles and 
anti-particles and hence not be positive definite. But for 
bosonic statistics it is non-negative. 

Note that, for a scalar field, the situation is completely 
different since in that case, the inner product is not pos- 
itive definite: (F*\F*) = - (F\F). 

V. BLACK HOLE ANALOGUE 

After introducing the notion of the effective geome- 
try we can now design an analogue of a black or white 
hole. To this end we move the medium with a constant 
velocity v in order to be able to correct (tune) the back- 
ground beam according to the resulting Dopplcr shift. 
Again the background solution, i.e., SI and 0, should be 
homogeneous if we want to avoid additional source terms 
for the linearized fields. (In the reference frame of the 
moving atoms, an inhomogencous background becomes 
time-dependent and thereby also causes a deviation from 
the dark state.) 

The only parameter left for influencing the effective 
geometry is the density p. In order to arrive at a station- 
ary effective metric, the density profile in the laboratory 
frame should be timc-indcpcndcnt p = p(x). In the rest 
frame of the fluid, this requirement implies a space-time 
dependence of p — p(x — vt). At a first glance, such 
a scenario seems to be inconsistent with a constant ve- 
locity v, but one could arrange a flow profile such as 
v = ve x — vyp'e y / p which, for a light beam at y = 0, re- 
produces these properties. As already mentioned at the 
end of Sec. IV A, an alternative possibility is to cause 
transitions \d) , \e) — > \a) , \b). 

Since we are still working with non-relativistic veloci- 
ties v -C 1, the rest frame of the medium and the labo- 
ratory are related by a Galilei transformation 

d d d d d 
Ot ot ox ox ox 

Having derived a covariant, i.e., coordinate-independent, 
representation of the effective action in Eq. (42), this 
transformation is completely equivalent to a correspond- 
ing change of the effective geometry 

7° - 7° , 7 1 - 7 1 +« 7 °. (58) 

The effective metric is then given by the well-known 
Painlcvc-Gullstrand-Lemaitre form [11] 

C = (l + 2K)(i , 2 _ 1/ ( 1 + 2 H))- ( 5Q ) 
The inverse metric simply reads 



<=( 1/(1 + r»-" 2 _l). <») 

Obviously, a horizon (g^ = 0) occurs for v 2 = 1/(1 + 
2H), which could be a relatively low velocity and perhaps 
experimentally accessible. 

A. Negative Effective Energy 

For stationary (in the laboratory frame) parameters 
H = H(x) and v = const one may construct a conserved 
energy (Noether theorem) of the beat fluctuations ip via 
Eq. (54). Since for a moving medium, the effective met- 
ric in Eq. (60) has off-diagonal elements, the resulting 
expression is more complicated than in Eq. (55). For 
the sake of convenience, we adopt the geometric-optics 
approximation u>, n ^> H'/H and obtain 

E = Jd X ^ U+ ^ + <T ' V ^^^ + VK U, (61) 

and, after diagonalization and normalization (ip\tp) = 1, 
the solutions for the effective energy assume the following 
form 

E eS = ^ (u + (l ± v VTT2N) {u + v «)) . (62) 

We observe that even the branch u> > of the disper- 
sion relation which corresponds to a positive energy in 
flat space-time can become negative beyond the horizon 
v > 1/Vl + 2H. This purely classical phenomenon - i.e., 
that the energy measured at infinity can become negative 
beyond the ergo-sphcrc g 00 = - occurs for real black 
holes as well and can be considered as the underlying 
reason for the mechanism of super-radiance, etc. 

Of course, the total energy of the system as derived by 
the total action in Eq. (1) is always positive. The modes 
with a negative effective energy (pseudo- or quasi-energy, 
cf. [12]) possess a total energy which is smaller than that 
of the background. In this regard a (classical) mixing of 
positive and negative (effective) energy modes is possible. 

B. Bogoliubov Coefficients 

If the effective metric possesses a horizon, one would 
expect the usual mixing of positive and negative energy 
solutions as governed by the Bogoliubov coefficients Qb 
and Pe defined via 

^ = a E ^T+pE^ E , (63) 

with the positive (tp E ) an d negative (i/>_e) energy 
modes, respectively, which are normalized 

(y£iv4 n ) = (^-eW-e) = 

(r E ut \r E ut ) = (r- E \r- E ) = i , (64) 
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and orthogonal 



)=0. 



(65) 



Owing the the positivity of the inner product in Sec. IV D 
the completeness relation has a plus sign instead of a 
minus as for the scalar field, i.e., 



\a E \ 2 + \P E \ 2 = 1. 



(66) 



Consequently we obtain the Fermi-Dirac factor for the 
scattering (Bogoliubov) coefficients (instead of the Bose- 
Einstcin distribution for scalar fields) 



\/3 E \ =e" E / K °\a E \ ~> \Pe\ 2 = 



(67) 



However, it should be emphasized here that this mode 
mixing is a priori a purely classical phenomenon and in- 
dependent of the quantum features (commutation rela- 
tions) - the fields i/> do not obey Fermi-Dirac statistics, 
see the next Section. Only if the quantum commutation 
relations assigned a physically reasonable particle inter- 
pretation to the modes ip E - as it is the case for a truly 
fermionic Dirac quantum field, for example, but not for 
the fields $± (see below) - one could infer the (quantum) 
Hawking radiation. 

The surface gravity of the effective horizon at v 2 = 
1/(1 + 2N) = c 2 low depends on the rate of change of the 
velocity of light in the laboratory frame across the hori- 
zon 



d(\v\ 



Cslow J 



dx 







horizon 


dx 



horizon 



an 

dx 



(68) 



horizon 



By comparison with Eq. (45), we observe that k s is of 
the same order of magnitude as the rest energy induced 
by the effective mass (remember the homogeneous Dirac 
equation (i-f°d t — ra)i\) = 0) 



TO 



1 



1 



an 



V1 + 2H 2 ^(1 + 2^)3 dx ' 



(69) 



As a result, the relevant mode-mixing effects - i.e., the 
Bogoliubov /3-coefficients - are not strongly suppressed 
by the effective mass. 



VI. COMMUTATION RELATIONS 

Having derived an effective metric which may exhibit a 
horizon, one is immediately led to the question of whether 
the system under consideration could be used to simulate 
the Hawking effect. As it will turn out, the answer is "no" 
- since the Hawking effect is a quantum effect, it is not 
sufficient to consider the wave equation, one also has to 



check the commutation relations which generate the zero- 
point fluctuations (the source of the Hawking radiation) 
according to the Heisenberg uncertainty principle. For 
convenience we shall transform back into the rest frame of 
the medium and assume a constant H for the calculations 
in this Section. 



A. Commutators 

Obviously the effective action derived above is intrin- 
sically different from the one of a charged scalar field, 
for example. To make the difference more explicit let 
us consider the effective (adiabatic limit) commutation 
relations following from Eq. (33). 

For any given time to, the equal-time commutation re- 
lations of the fields <l± vanish. Since the equations of 
motion do not mix $± with <&±, this remains true for all 
times 



#±(t,*),$±(f,a;') - & ± (t,x),& ± (t',x') =0. (70) 



According to Eq. (33) the canonical conjugated momenta 
are i&+ and i[l + 2H]$1, respectively, and hence we ob- 
tain 



$ + (t,x),& + (t,x')} =S(x-x'), 



and 



6(x-x') 



1 + 2H 

The remaining (equal-time) commutators vanish 
& + (t,x),$-(t,x')] = [$+(*,*), $*-(*, a/) 



(71) 



(72) 



0, (73) 



and the commutation relations for the time-derivatives 
of the fields can be inferred from the equations of motion 
<!>+ + = and (1 + 2K)$_ + $' + = 0. 

Remembering the definition of the effective two- 
component spinor in Eq. (40) the above relations can 
be cast into the compact form 

t> A {t,x),i> B {t',x')\ - W A (t,x),^ B (t',x')] =0, (74) 



as well as 



^ A (t,x),i} B (t,x') =5 A b5(x-x'). 



(75) 



Since the beat fluctuation of the electromagnetic field 
(coupled to the medium) do not obey the Pauli exclusion 
principle, one cannot fill the Dirac sea consisting of all 
negative (effective) energy (in flat space-time) states and 
thereby define a new vacuum state - as it is possible for 
fermionic quantum fields. 
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B. Comparison with Other Fields 



C. Particle Creation 



Let us compare the above commutation relations with 
those of a (1+1 dimensional) Schrodinger field ip 



as well as 



^{t.x),^^^ 1 )] = 5{x - x 1 ) , 



= 0, (76) 



(77) 



on the one hand and and with the commutators of a (1+1 
dimensional) charged scalar field <f) 



j>(t,x),${t',x')\ = [ft(t,x),ft(t',x') 
j>{t,x),ft{t,x')] = \ft{t,x),${t,x') 



0, 
0, 



as well as 



<f>(t, x), dt$(t, x') = i5(x — x') , 



(78) 



(79) 



on the other hand. In the latter case (charged scalar 
field 0), the equation of motion can mix positive and 
negative frequencies and thereby lead to particle produc- 
tion - whereas in the former situation (Schrodinger field 
■0), the number of particles is conserved. This difference 
becomes more evident when one decomposes the fields 
into real (self-adjoint) and imaginary (anti-self-adjoint) 
parts. For tp, the independent canonical conjugated vari- 
ables are $iip an d 3V> ~ whereas for cf), they are and 
(as well as 30 and 30). 

Obviously, the commutation relations of the fields <fr± 
are clearly inconsistent with those of a charged scalar field 
and show more similarity to the (bosonic) Schrodinger 
field. Therefore, the system under consideration cannot 
serve as a true analog for the quantum effects in the pres- 
ence of a black hole horizon - such as Hawking radiation 
- although it reproduces all classical phenomena. 

Since the fields $± describe fluctuations of the electro- 
magnetic field, it is also clear that they do not obey the 
fermionic (anti) commutation rules 

{j> A (t,x),j> B (t,x')} = {j> A (t,x),j> B (t,x')} =0, (80) 

as well as 

{j> A (t,x),ipUt,x')} =S AB S(x-x'). (81) 

An effective Dirac field satisfying bosonic commutation 
relations might seem rather strange in view of the spin- 
statistics theorem. Indeed, one key ingredient needed 
for the derivation of this theorem, the spectral condition 
(which is one of the Wightman axioms) , is not satisfied in 
our case, since the effective energy can become negative 
owing to the huge total energy of the background field, 
see also Sec. IV C. 



In order to answer the question of whether there is 
any particle creation at all in the described slow-light 
system, one has to clarify the notion of (quasi)particlcs 
to be created (or not) and to specify the corresponding 
(in/out) vacuum state. 

For example, an appropriate initial state |in), which 
is a coherent state in terms of the fundamental creation 
and annihilation operators of the electromagnetic field, 
could be chosen such that it is annihilated by all fields 
*±, 



Vx : $ + (t in ,x) |in) = $_(t in ,a;) |in) = 0. 



(82) 



This is possible because the fields <I>± are purely decom- 
posed of positive frequency parts of the electromagnetic 
field, i.e., the annihilators, cf. Eq. (7). If the effective 
Hamiltonian of the fields $± (in an asymptotically flat 
region, i.e., for a homogeneous medium at rest) is given 
by a non-negative bilinear form such as 



T>i> 



(83) 



<3> + | in) = $_ | in) = is indeed the (or at least one) 



with T> denoting a (differential) operator, the state 
|in) 

ground stated 

In this case the initial (vacuum) state is annihilated by 
the fields <l± at all times 



\/t,x : <f>+{t,x) |in) = $-{t,x) |in) = 0, 



(84) 



as the time-evolution does not mix $± with and 
there is no particle creation. 

For another initial (vacuum) state (e.g., a squeezed 
state) and a different particle concept, 



/ 



0. 



(85) 



however, some effects of (quasi)particle creation might 
occur. These phenomena could be tested by sending in a 
(multi-mode) squeezed state and comparing the number 
of photons per mode in the in- and out-states. 

Another possible source for (quasi) particle creation is 
the finite life-time of the atomic state |c) as represented 
by the effective decay rate T. Realistically, this decay 
corresponds to some spontaneous emission process gen- 
erated by the quantum fluctuations of the electromag- 
netic field, for example. Consistent with the fluctuation- 
dissipation theorem this coupling also introduces (quan- 
tum) noise, which is not included in our treatment and 
could possibly lead to particle creation. However, this is 
clearly a pure trans-Planckian effect and cannot be in- 
terpreted as Hawking radiation. 



''Therefore, it cannot be the equivalent of the Israel-Hartle- 
Hawking [13] state, in which the Hawking radiation is some- 
what hidden by the fact that there is no net energy flux. 
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VII. DISPERSION RELATION 

Although slow light cannot be used to simulate the 
Hawking effect it can reproduce various classical effects 
associated to horizons**, such as mode mixing and the 
associated Bogoliubov coefficients, see Sec. VB. In view 
of the red- or blue-shift near the horizon deviations from 
the linear dispersion relation have to be taken into ac- 
count, cf. [16]. With the choice in Eq. (23) the dispersion 
relation in Sec. Ill A simplifies because of X(u) = Y(uj), 
and we obtain for a medium at rest, cf. Figs. 1 and 2 




-10 -5 5 



FIG. 1. One branch of the dispersion relation of the (Ffield 
in Eq. (86). Frequency lo and wave-number k are plotted in 
units of the Rabi frequency Q for K = 10 and T/Q = 1/10. 
These values (of order one) are but illustrative and chosen 
in order to resolve the characteristic features in one figure - 
realistically the orders of magnitude arc different. The imag- 
inary part describes the absorption and does not change sig- 
nificantly in the limit r J, 0. For very large as well as for 
very small lo the medium becomes transparent. The steep 
slope within the transparency window u> <C corresponds to 
the reduced propagation velocity - whereas the effect of the 
medium for large lo is negligible. As one can observe, the 
anomalous frequency solutions lo > Q are separated from the 
normal ones lo < Q by a large region of absorption. 

We observe two major differences between the disper- 
sion relation above and that for the sonic black hole 



analogs, for example in Bosc-Einstcin condensates (see 
[17] and Sec. IX A) with 

, 2 = Cs 2 onn /(l + (V) , (87) 

where £ denotes the so-called healing length and provides 
a wave-number cut-off, cf. Fig. 3. 



10 - 



5 - 



3 o - 



-5 - 
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K 

FIG. 2. The real part of the dispersion relation in Fig. 1 as 
lo vs. k with the same values. One can easily recognize that 
the first deviation from the linear dispersion relation at lo <C 
is "subluminal" - although it becomes finally "superluminal" 
for lo S> fi. The solutions with an anomalous (negative or 
even infinite) group velocity lie completely in the absorptive 
region, cf. Fig. 1. 

The sonic black hole analogs generate a deviation from 
the linear dispersion relation via the spatial dependence 
(k) and, consequently, for each value of the wave-number 
k there exist two possible solutions for the frequency (±ui 
for a medium at rest). In contrast, for the black hole 
analogs based on slow light the deviation is mainly^ 
caused by the (non-local) temporal dependence. (This 
remains true for all dielectric/optical black hole analogs, 
cf. [6,8].) As a result, one has two values of k for each 
value of u), but can have more than two solutions for u> 
for some values of n. Even though these anomalous solu- 
tions for lo are separated from the normal ones by a rela- 
tively large region of absorption, it would be interesting 
to see under which circumstances this peculiar behavior 
may give raise to additional effects (such as mode mixing, 
etc.). 



^Of course, the finite interatomic distance results in a devi- 
"Other systems which are potentially capable of simulat- atkm from thg Unear dispersion relation too> but the cut _ off 

ing those classical effects with present-day technology are dis- given by thfi RaW frequency is usually reached earlier . 

cussed in Refs. [14] and [15]. 
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Another major difference between the dispersion re- 
lations (86) and (87) is that the sonic dispersion rela- 
tion (87) is "superluminal" /supersonic for large wave- 
numbers Wgroup = duj/dn > c soun d for £k ^t- 1 whereas 
the slow-light dispersion relation (86) is "subluminal" 
^group = dio/dn < l/vl + 2H within the transparency 
window, say \u\ < 0/2, but |w| ft. For very large 
frequencies one recovers the speed of light in vac- 

uum u> = k — although this limit is totally outside the 
region of applicability of our approximations. 




-2-1 1 

K 



FIG. 3. One branch of the dispersion relation of (zero) 
sound waves in Bose-Einstein condensates at rest, cf. Eq. (87), 
in arbitrary units. If the condensate is moving the vari- 
ous K-solutions for a given frequency ui in the laboratory 
frame can be found by the points of intersection with straight 
lines as determined by Eq. (88). For a subsonic velocity 
v < c soun d, there is only one solution, denoted by s+, which 
has a small wave-number and a positive pseudo-norm, i.e., a 
positive c.j nuid / srcst _ framc (assuming t^iab-framc > 0). For su- 
personic velocities, on the other hand, i.e., beyond the hori- 
zon, there are three possible solutions - one with a small 
wave-number and a negative pseudo-norm (s-) as well as two 
others with large wave-numbers and positive (1+) and nega- 
tive (1-) pseudo-norm, respectively. The mixing between these 
modes at the horizon generates the Hawking radiation (s+). 

VIII. PROBLEMS OF SLOW LIGHT 

The direct (naive) way to use the most common set-up 
for slow- light experiments - i.e., a strong control beam 
and a weak (perpendicular) probe beam - in order to 
build a black hole analog goes along with a number of 
(somewhat related) difficulties listed below. Whereas the 
first three obstacles are can be avoided by the arrange- 
ment proposed in this article, the fourth one persists - 
indicating that this system is a classical, but not a quan- 



tum analogue of a black hole. 




-40 -20 20 



FIG. 4. One branch of the dispersion relation of a slow-light 
pulse (in the usual set-up) k 2 = lu 2 [1 + p(u + lu ) — p(u) — lu )] 
where p(uo) = 2K(fi 2 / 'wo)w / '(to 2 - O, 2 + iTui), see e.g., [1,10], 
in units of the Rabi frequency O for ujo/Q = 20, T/Q = 1/2, 
and K = 5. Again, these unrealistic values have been cho- 
sen in order to illustrate the chracteristic features. For more 
realistic values the peaks would be more pronounced, the 
transpacency windows \u) ± u>o| <fC f2 narrower, and the slope 
inside them steeper, etc., but the main structure remains. For 
|w±u;o| S> the influence of the medium is negligible. Within 
the transpacency windows \u ± u>o| -C fi, the steep slope in- 
dicates a reduced group velocity and the solutions with an 
anomalous group velocity \uj ± u>o| = 0(Q) lie inside the ab- 
sorptive regions. 



A. Frequency Window 

Light pulses (of the probe beam) are only slowed down 
drastically - or may propagate at all - in an extremely 
narrow frequency window in the optical or near-optical 
regime. But the frequency of the particles constituting 
the Hawking radiation cannot be much larger than the 
surface gravity (e.g., the gradient of the fluid's velocity) 
which makes an experimental verification in this way very 
unlikely. 



B. Doppler Shift 

In a stationary medium, the frequency as measured in 
the laboratory frame is conserved - but the frequency 
in the atom's rest frame changes as soon as the velocity 
of the medium (Doppler shift) or the wave-number (red- 
shift) varies (which necessarily happens near the hori- 
zon). Hence the beam will leave the narrow frequency 
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window - which is generated by the (moving) atoms - in 
general. 

C. Group and Phase Velocity 

Since the group and the phase velocity of the probe 
beam are extremely different w gr0 up -C v p hasc w 1, it is 
not possible to describe its dynamics by an effective local 
wave equation resembling a scalar field in a curved space- 
time. 




-40 -20 20 
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FIG. 5. The real part of the dispersion relation in Fig. 4 
as lo vs. k with the same values. The additional line demon- 
strates the slope corresponding to a motion of the medium 
with the reduced group velocity as in Fig. 3. Obviously, there 
can be no mixing of positive and negative pseudo-norms via 
the usual mechanism sketched below Fig. 3 in this case. Even 
though the peaks can be much higher for small F and thereby 
could possibly intersect with the straight line, the resulting 
solutions would lie completely in the region of strong absorb- 
tion (cf. Fig. 4) and therfore do certainly not model Hawking 
radiation. 



D. Positive/Negative Frequency-Mixing 

In order to obtain particle creation, one has to have a 
mixing of positive and negative frequencies, or, more ac- 
curately, positive and negative pseudo-norm (as induced 
by the inner product, cf. Sec. VB) solutions. In a station- 
ary flowing medium (as used for the black hole analogs), 
this can occur by tilting the dispersion relation due to 
the Doppler effect caused by the velocity of the medium 

<^lab — frame ^fluid's rest — frame "f" ^medium ^ • (88) 

As soon as the velocity on the medium exceeds \co/k\, 
i.e., the phase velocity, a mixing of positive and negative 



frequencies (in the fluid's rest frame) becomes possible, 
cf. [9]. However, since the phase velocity of the slow-light 
pulse is basically the same as in vacuum, this mechanism 
does not work in this situation and, consequently, there 
is no particle creation. 

IX. COMPARISON WITH OTHER SYSTEMS 

One of the main points of the present article is the 
observation that an appropriate wave equation and the 
resulting effective geometry of a black hole analog is not 
enough for predicting Hawking radiation. Although all 
the classical effects can be reproduced in such a situation, 
the adequate simulation of the quantum effects requires 
the correct commutation relations as well. 

In view of this observation one might wonder whether 
this is actually the case for the currently discussed (e.g., 
sonic/acoustic and dielectric/optical) black hole analogs. 
In the following we shall deal with this question for two 
representative examples, for which the commutation re- 
lations can be derived easily. 

A. Bose-Einstein Condensates 

The dynamics of Bose-Einstein condensates are to 
a very good approximation described by the Gross- 
Pitaevskii equation 

^=[-^ + v ( r )+V 2 H\ 2 ^, (89) 

where ip denotes the mean-held amplitude, m the mass 
of the bosons, V an external (trapping) potential, and rj 
is the scattering parameter governing the two-body re- 
pulsion of the constituents. Inserting the eikonal ansatz 
(Madelung representation) , 

^ = ^e* s , (90) 

and introducing the (mean-held) velocity v = VS/m, 
one obtains the equation of continuity g + V(gw) and 
the equivalent of the Bernoulli or the Hamilton- Jacobi 
equation 

> + r + *^-^- «»> 

Within the Thomas-Fermi approximation, one neglects 
the quantum potential, i.e., the term on the l.h.s., and 
hence recovers the usual equations of fluid dynamics, see 
also [17]. The linearization around a given (stationary) 
background profile qq and So — * vq yields the well-known 
wave equation 

{d t + V • v ) (d t + v • V) SS = V- Vg VSS . (92) 
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The commutation relations of SS, which we are interested 
in, can be derived from the commutator of the fundamen- 
tal fields 



${t,r),ft(t,r') = S 3 (r-r'). 



(93) 



Inserting the linearization of ip = \fgcxp(iS) around a 
classical background via g = gq + Sg and S = So + SS we 
obtain (note that g = g^ and S = S^) 



Sg{t,r),SS(t,r') = i5 3 (r-r'). 



(94) 



The relation between Sg and SS follows from Eq. (91) in 
the Thomas-Fermi approximation 



5g= — -(d t +v -V)SS. 



(95) 



Hence Sg is indeed the (negative) canonical conjugated 
momentum to SS - provided that one inserts the constant 
factor rj 2 correctly into the (effective) action - and the 
commutation relations are equivalent (within the used 
approximation) to those of a quantum field in a curved 
space-time. 



B. Non-Dispersive Dielectric Media 

As another example we study non-dispersive and linear 
dielectric media, see e.g. [6]. For a medium at rest the 
fundamental Lagrangian describing the electromagnetic 
field, the dynamics of the medium (£[P]), as well as their 
mutual interaction (E • P) is given by 

C = ^(E 2 -B 2 )+E-P + C[P}. (96) 

Accordingly, using the temporal gauge and introducing 
the vector potential via E — d t A and B = — V x A, the 
canonical momentum is just the electric displacement 



n = d = e + p. 



(97) 



Performing basically the same steps as in Sec. IV A we 
may integrate out the degrees of freedom associated to 
the medium P and thereby arrive at an effective (low- 
energy) action for the (macroscopic) electromagnetic field 
alone, cf. [6]. But, in contrast to the highly resonant 
behavior of P in slow-light systems, non-dispersive media 
respond adiabatically with a constant susceptibility \ = 
e — 1, i.e., P = \E and thus II = D = eE, to the 
external field (at sufficiently low frequencies), cf. [6]. 

If the (non-dispersive) medium is moving with the ve- 
locity f3 the electric and magnetic fields get mixed and 
one obtains 



n = D = eE + (e - l)B x /3 + C(/3 2 



(98) 



Again, the commutation relations fit to an effective- 
metric description - which is not completely surprising 
because the effective action has the same form as in 
curved space-times, cf. [6]. 



X. DISCUSSION 

Let us summarize: The naive application of slow 
light (i.e., the most common set-up) in order to create 
a black hole analog goes along with several problems, 
cf. Sec. VIII. With the scenario proposed in this article, 
the problems associated to the classical wave equation 
can be solved and it is - at least in principle - possible 
to create a (classical) black hole analog for the <& field. 
At low wave-number, the corresponding dispersion rela- 
tion represents a quadratic relation between k and ui, and 
can thus be written in terms of an effective metric. If the 
fluid is in motion, this low wave-number equation can be 
changed into a black hole type wave equation. 

However, this classical black hole analog does not re- 
produce the expected quantum effects - such as Hawking 
radiation**. In order to simulate the Hawking effect, it is 
not sufficient to design a system with an equivalent effec- 
tive equation of motion - the commutation relations have 
to match as well. This is indeed the case for the sonic 
black hole analogs in Bose-Einstein condensates and non- 
dispersive dielectric black hole analogs - but for sound 
waves in more complicated systems, for example, it is 
not immediately obvious. 

Nevertheless, in the scenario described in this article, 
the field <& governing the beat fluctuations of an electro- 
magnetic background field obeys the same equation of 
motion as in the presence of a horizon and hence can be 
used to model several classical effects associated to black 
holes - for example the mode mixing at the horizon as 
described by the Bogoliubov coefficients, see Sec. VB. 
One way of measuring the Bogoliubov coefficients could 
be to send in a "classical" pulse above the background 
- i.e., a particular coherent sate in terms of the funda- 
mental electromagnetic field - and compare it with the 
outcoming pulse. As another (more fancy) possibility one 
might think of a multi-mode squeezed state - which in 
some sense simulates the vacuum fluctuations which are 
transformed into quasi-particles by the mode mixing. 

However, one should bear in mind that, as the 
wave-packets propagate away from the horizon and get 
strongly blue-shifted, they eventually reach the regime 
where the concept of the effective geometry breaks down 
and effects like dispersion, non-locality (in time) of the 
effective action, and, finally, absorption become rele- 
vant. For a reasonably clean interpretation, therefore, 
one should investigate the scattering of the wave-packets 
not too far away from the horizon. 



H This conclusion applies in the same way to the scenario 
proposed in Ref. [7], where the Schwarzschild metric is sim- 
ulated by a medium at rest with the horizon corresponding 
to a singularity in the effective refractive index. Such static 
analogs of the Schwarzschild geometry (see also [18]) go along 
with further problems [19]. 
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A. Miles Instability 

Another interesting classical effect is related to the neg- 
ative parts of the energy in Eq. (62). Since a conserved 
positive definite energy functional of the linearized per- 
turbations would demonstrate linear stability, the neg- 
ative contribution in Eq. (62) can be interpreted as an 
indicator for a potential instability (e.g., super-radiance) 
- provided a suitable coupling between positive and neg- 
ative (effective) energy modes. 

As an example, let us assume that the "superlumi- 
nally" flowing v > l/vl + 2H slow-light medium inter- 
acts with the environment in the laboratory frame via a 
friction term such as Td t Q (with possible spatial deriva- 
tives). For small ui and k the resulting dissipation alters 
the dispersion relation via 

^ + v K ) 2 = y ^-zluT( K ) 1 (99) 

with the potentially K-dependent (additional spatial 
derivatives) friction term T(k) describing the interaction 
of the $-ficld with the environment at rest. 

For small T the imaginary part of the solutions for the 
frequency uj (assuming a real wave-number k € R) reads 

s>(w) = (i ± vVT+m) . (ioo) 

Consequently, beyond the horizon v > l/vT+~2K one 
of the allowed frequency solutions acquires a positive 
imaginary part and thus the dissipation (interaction with 
the environment) generates an instability. Note that the 
relative velocity v > \j\J\ + 2H between the slow-light 
medium and the environment (at rest) is crucial since a 
friction term like T(dt + vd x )& — -> + vk)T would of 
course not lead to any instability. 

This instability is somewhat analogous to the Miles in- 
stability [20] generating surface waves in water by wind 
blowing over it. In Rcf. [15], this phenomenon is called 
thermodynamic instability since it occurs when the free 
energy of the medium acquires negative parts in the 
frame of the environment. 



B. Outlook 

Apart from the aforementioned experiments there are 
many more conceivable tests one could perform with the 
proposed classical black hole analog based on slow light. 
A more drastic way of investigating the interior struc- 
ture of the sample (than the mere comparison of the 
in- and out-states) could be to freeze the dark state by 
completely switching off the background field and take 
a "snap-shot" of the state of the atoms by illuminat- 
ing them with strong laser beams with frequencies cor- 
responding to certain atomic transitions and measuring 
the absorption. 



Furthermore it would be interesting to investigate the 
influence of the anomalous frequency solutions of the dis- 
persion relation generated by the non-local temporal de- 
pendence (cf. Sec. VII), for example, on additional mode- 
mixing. This question is relevant for more general (non- 
dispersive) dielectric black hole analogs and might also 
lead to some insight into the trans-Planckian problem. 
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